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Abstract 

Using numerical simulations of quenched SU(2) gauge theory we demonstrate 
that an external magnetic field leads to spontaneous generation of quark 
condensates with quantum numbers of electrically charged p mesons if the 
strength of the magnetic field exceeds the critical value eB c = 0.927(77) GeV 2 
or B c = (1.56 ± 0.13) • 10 16 Tesla. The condensation of the charged p mesons 
in strong magnetic field is a key feature of the magnetic-field-induced elec- 
tromagnetic superconductivity of the vacuum. 

Keywords: Quantum Chromodynamics, Strong Magnetic Field, Phase 
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Effects caused by very strong magnetic fields attract increasing inter- 
est motivated by the fact that the hadron-scale strong magnetic fields may 
emerge in the heavy-ion collisions at the Relativistic Heavy-Ion Collider 
at Brookhaven National Laboratory and at the Large Hadron Collider at 
CERN [T] . Such fields may, presumably, have arisen in the early Universe [2] . 
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The strong magnetic field causes exotic effects in hot quark-gluon matter, 
the well-known example is the chiral magnetic effect [3J. In the absence of 
matter the magnetic field background leads also to unusual effects like the 
magnetic catalysis [4J , shift of finite-temperature transitions in QCD [5] and 
anisotropic conductivity [5]. 

It was recently suggested that a sufficiently strong external magnetic field 
should turn the vacuum into an electromagnetic superconductor (TJ |8]. The 
superconductivity emerges due to spontaneous condensation of electrically 
charged vector particles, p mesons, if the magnetic field exceeds the critical 
strength 

eB c = m 2 p ^ 0.6 GeV 2 , B c = m 2 p /e « 10 16 T, (1) 

where m p = 775.5 MeV is the mass of the p meson and e is the elementary 
electric charge. In terms of the up-quark (u) and down-quark (d) fields the 
suggested p- meson condensate should have the following form [Hj: 

(ujid) = -i(uj 2 d) = p , (uj 3 d) = (uj d) = , (2) 

where p = p(x L ) is a certain complex-valued periodic function of the coor- 
dinates x L = (x 1 , x 2 ) of a plane which is perpendicular to the magnetic field 
B = (0,0,B). 

The superconducting vacuum should have many unusual features. Firstly, 
no matter is required to create the superconductor so that the electromag- 
netic superconductivity appears literally "from nothing". Secondly, the su- 
perconductivity is anisotropic so that the vacuum acts as a superconductor 
along the magnetic-field axis only. Thirdly, the superconductivity is inhomo- 
geneous because the p-meson condensate is not uniform in the 5-transverse 
directions due to the presence a new type of topological defects, the p vor- 
tices. Fourthly, the net electric charge of the superconducting vacuum is zero 
despite of the presence of the charged condensates (|2| [S] . 

The spontaneous generation of the p-meson condensate ^ - which plays 
a role of the Cooper pair condensate in the conventional superconductivity 
- is the key feature of the vacuum superconductor mechanism [H [8] . 

The appearance of the p-meson condensate was found analytically in a 
phenomenological model based on the vector meson dominance (VMD) [7j, in 
the Nambu-Jona-Lasinio (NJL) model jH], and in holographical approaches 
based on gauge/gravity duality [9]. We use numerical simulations of the 
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lattice gauge theory to demonstrate that strong magnetic field indeed leads 
to emergence of the superconducting condensate of the charged p mesons. 

In QCD the charged p meson field is identified with p M = u^^d. The 
condensation pattern ^ corresponds to the condensate of the p mesons 
with the spins aligned along the axis of the magnetic field, with the s z = +1 
projection of the spin onto the z axis. It is convenient to introduce two 
combinations of the negatively-charged p-meson fieldfl 

P± = g (Pi ± if*) = 2^(Ti ± l ^) d , (3) 

which correspond to the spin projections s z = ±1, respectively. Indeed, 
according to the simplified arguments of Ref. [7j, the invariant masses M n>Sz 
of the p mesons states in the magnetic field B should behave as follows, 

Ml Sz = ml + (l + 2n-2s z )\eB\, (4) 

with the nonnegative integer n and the spin projection onto the z axis s z = 
— 1,0, +1. The ground state is identified with the quantum numbers n = 
and s z = +1, and the charged p mesons should get condensed, Mq +1 < 0, if 
the magnetic field exceeds the critical value ([T]). 

The simplest way to check numerically the possible appearance of the 
superconducting condensate (p| is to calculate the equal-time correlation 
function along the direction of the magnetic field, 

G±{z) = (pl(0) P ±(z)) , (5) 

where the separation in the transverse coordinates of the two probes is set 
to zero, x 1 - = (0,0). The long-distance behavior of the s z = +1 correlation 
function ^ should expose the expected emergence of the condensate §2§ due 
to the factorization property: 

lim N _ 00 G + (^) = |(p)| 2 , (6) 



2 One can equivalently work with the positively-charged /9-meson fields, <i(7i ± 172)^. 
The magnitudes of the positive and negative condensate are equivalent because the vacuum 
state is electrically neutral. Our results on vacuum condensation are the same for positive 
and negatively charged operators. 
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while the p mesons with the opposite orientation of the spins, s z = — 1, 
should not be condensectB 



lim 



G-{z) = 0. 



(7) 



We calculate the correlation functions ([5]) numerically, using lattice Monte- 
Carlo simulations of quenched SU(2) lattice gauge theory following numerical 
setup of Ref . [6] • The quark fields are introduced by the overlap lattice Dirac 
operator T> with exact chiral symmetry [10] . The correlation function ([5| is 
a linear combination of the current- current correlators in the vector meson 
channel. The vector correlator is represented in terms of Dirac propagators 
in fixed background of Abelian and non-Abelian gauge fields and is then 
averaged over an equilibrium ensemble of non-Abelian gauge fields 



where Sym [AJ\ is the lattice action for gluons A^. 

A uniform time-independent magnetic field B is introduced into the Dirac 
operator Vj for the flavor / = u, d in a standard way by substituting the 
su (2)- valued vector potential A^ with the w(2)-valued potential A^ij — > 
A^ij + S^qf F^v x u /2, where q$ is the electric charge of the corresponding 
quark, q u = +2e/3, qd = — e/3, and i,j are color indices. We also introduce 
an additional twist for fermions in order to account for periodic boundary 
conditions in spatial directions [TTI [T2"] . For technical reasons the bare quark 
mass mo is fixed at a small value am® = 0.01, where a is the lattice spacing. 
The vector correlation functions depend very weakly on the bare quark mass 
if they are calculated with the help of the overlap Dirac operator [T3] . 

Our numerical approach is done in two complimentary ways. Firstly, 
we study in details the superconducting condensate for eleven values of the 
magnetic field B at a relatively small symmetric 14 4 lattice at a single lattice 
gauge coupling (3 = 3.281. These parameters correspond to the physical vol- 
ume of the lattice is L 4 = (1.44 fm) 4 and the lattice spacing a = 0.103 fm [14] . 



3 The condensed component is identified by the positiveness of the product eBs z . If the 
magnetic field is reversed and becomes negative, eB < 0, then the condensed component 
corresponds to s z — — 1 while the component with s z — +1 is not condensed. 



( u (x) 7^ (x) d (y) ^ v u (y) ) 
■ f DA„e- SYM[A ^Tr ( — 




(8) 
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Then we use an heuristic fitting method to find the superconducting conden- 
sate without taking a long-range limit ^ because the factorization ^ does 
not work well in too small volume. Secondly, we consider a set of lattices of 
various physical volumes L 4 and four values of magnetic field strengths B, 
and then utilize a conventional fitting procedure to extract the condensate 
77 = 1](L). The extrapolation to the infinite volume, L = la — > 00, shows that 
these two methods give the same results. In both approaches the ultraviolet 
lattice artifacts are reduced with the help of the tadpole-improved Symanzik 
gluon action [13] . 

In our first approach we use 30 configurations of the gluon gauge field for 
each value of the background magnetic field. The periodicity of the lattice 
leads to the quantization (k 6 Z) of the magnetic field, 

3 . 2n 

B = kB minj eB min = — = 0.354 GeV 2 , (9) 
because of the requirement 

J d 2 x ± q f B E Z for f = u,d. (10) 

In Eq. ([9]) the integer k = 0, 1, ... , L 2 s /2 determines the number of elementary 
magnetic fluxes which pass through the boundary of the lattice in the (x l , x 2 ) 
plane. 

The maximal possible value of the fluxes k = l 2 /2 = 98 corresponds to 
an extremely large magnetic field with the magnetic length Lb ~ (eB)^ 1 ^ 2 
being of the order of the lattice spacing, Lb ~ a- In order to avoid associated 
ultraviolet artifacts, in our simulations we limit the maximal value of the 
fluxes by k max = 10 I 2 / '2, so that our maximal magnetic field, el? max = 
3.54 GeV 2 is much larger than the expected critical magnetic field ([!]). 

In Fig. [l] we show correlator ^ for a few relatively small values of the 
magnetic field. As we have anticipated, the s z = ± 1 correlators split in 
the external magnetic field [in Fig. [I] we show both spin orientations for 
B = -B m in]. The observed splitting reflects the change in the relevant lowest 
energies (El), M± = m 2 =1= \eB\, so that the expected hierarchy of the masses, 
M + <m p <M_, is encoded in the slopes of the correlators 

The splitting of the s z = ±1 masses was also found in SU (3) lattice gauge the- 
ory at weak magnetic fields [33]. We have checked that the no- condensation 
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Figure 1: The G+ correlator (s z = + 1) at eB=Q, 0.35, 0.71 GeV 2 and the G_ correlator 
(s z = — 1) at eB=0.35 GeV 2 in the insulator (weak magnetic field) phase (in lattice units). 



The lines illustrate the best fits by the heuristic function (111 



property ^ for the s z = — 1 correlator G- is valid for all studied values of 
the magnetic field. 

In the weak field region, B < B c , the long distance behavior of the corre- 
lator G + is expected to be proportional to the function e~ M ' z ', or cosh[(|^| — 
L/2)fi] in a finite volume with periodic boundary conditions. Here fi is a 
massive parameter and L = la is the physical lattice size. 

We have found, however, that our numerical data for the correlators ^ 
are consistent with the cosh-like behavior only in a very narrow interval of 
the coordinate z. Therefore, in at B < B c we used a heuristic fit function 

G^ weak = C weak e-^ L / 2 coshXN - L/2)\ , (11) 

which reduces to the exponential ( "no-condensate" ) behavior in the thermo- 



dynamic limit L — > oo. In Eq. (11) [i > 0, 7 > and C wea k > are the 
fitting parameters. 



The function (11) which works surprisingly well at weak values of the 
magnetic field ^ with k = 0,1,2. The values of ^/d.o.f. are shown in 
Table [TJ The heuristic function fits nicely our numerical data for the G + 
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2.1 


1.9 


1.1 


0.23 


0.25 


0.37 



Table 1: Values x 2 /d.o.f for the fit functions in the regions of weak (11) and strong (12) 
magnetic field (separated by the double vertical line) vs the magnetic field eB for 14r 
lattice. 



correlator (|5j in the weak field domain while at larger fields (B ^ 1.06 GeV 2 ) 
the fit gives unacceptably high values of x 2 /d-o.f. The best fit values for 
the parameter 7 are 7 ~ 5 ... 7. The k — 0, 1 (k=2) fits exclude the short 
distance separations z ^ a (z ^ 2a) and their periodic mirrors. The data 
and the best fits are shown in Fig. [TJ 

We have found that in the high-strength region, eB > 1 GeV 2 , the numer- 
ical data for the correlator G + (z) can well be described by another heuristic 
function 

£,fit,strong/ 



V{z) = C s tron g e"^ /2 cosh[/i(| 2 | - L/2)] 



(12) 
(13) 

a (and its 



for all separations z excluding the ultraviolet region with z ^ 
periodic mirror). In Eqs. (12) and (13) p > 0, r\ > and C n i g h < are the 
fitting parameters. The G + correlator in the strong magnetic field region 
and the corresponding best fits (12) are shown in Fig. |2l 



The parameter rj in the fitting function (12) plays a role of the charged 
p-meson condensate, 77 = \(p)\, because in the thermodynamic limit, L — > 00, 
the function V(z) reduces to an exponential e - ^' so that 

-(fit, strong / 



lim lim G H 

Z— KX) l s — >oo 



\z) = rf 



(14) 



Thus, the fits of the G + correlators provide us with the values of the conden- 
sate of the charged p mesons, Fig. [3} The corresponding values of x 2 /d.o.f. 
are shown in Table [TJ In the weak field domain the fitting does not converge 
properly due to presence of flat directions in the fitting parameter space. 

In the weak field region the p meson condensate rj = \(p)\ vanishes, while 
at higher values of the magnetic field the condensate deviates spontaneously 
from zero signaling the presence of the superconducting phase. We find that 
the dependence of the p-meson condensate on the magnetic field can be 
described by the linear function: 



77(5) = C p ■ {eB - eB c 



B>B r . 



(15) 
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Figure 2: The G+ correlator ^ in the superconducting (strong magnetic field) phase. 
The lines are the best fits by Eq. (fl2|). 



The best linear fit (shown by the dashed line in Fig. [3]) of the condensate al- 
lows us to determine the critical magnetic field of the insulator-superconductor 
transition, 

eB c = 0.924(77) GeV 2 , (16) 

or B c = (1.56 ± 0.13) ■ 10 16 T, in satisfactory agreement with the theoretical 
relation eB c = m 2 p , for the quenched mass of the p meson in SU(2) 



lattice gauge theory, m p ~ 1.1 GeV [17] . The prefactor in Eq. (15) is C p = 
7.5(5) MeV. Notice that in our quenched model the exponent in Eq. (15) 
is v = 1 while the mean field methods both in the bosonic VMD model [7] 
and in the fermionic NJL model [8] predict v — 1/2 (so that theoretically 
rj ~ - B c for B ^ B c ). 

The behaviour of the fitting parameter /i in the fitting functions at both 



sides of the critical phase transition (11) and (12), (13), are shown in Fig. 4 



The unusual forms of the fit functions (11 ) and (12) is used to absorb the 
finite volume effects. In order to support our small-volume results we have 
performed an infinite-volume extrapolation of the condensate r] obtained at 
larger lattices (/ = 17. . . 21 with L w 1.65 . . . 2.2 fm). The condensate was 
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Figure 3: The superconducting condensate r\ — of the charged p mesons as the 

function of the magnetic field B. The green points correspond to the condensate calculated 
for small lattice 14 4 , while the blue squares represent the data extrapolated to an infinite 



volume L — > oo. The dashed blue line is the fit by the linear function ( 15 ). The red arrow 



marks the point of the insulator-superconductor phase transition (j 1 €3 ) 
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Figure 4: The massive parameter /i corresponding to the best fits (11) and (12), (13) 
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Table 2: The parameters of the lattices used in the thermodynamic extrapolation for 
nonzero values of B and the corresponding best fit parameters obtained with the help of 
Eq. (17). The data is visualised in Fig.pSl 



obtained by fitting of the numerical data by the standard function, 

G%{z,L) = Acosh[m{z - l s /2)] + V 2 {L) . (17) 

where A, m and 77 are the fitting parameters. The fitting parameters and 
parameters of the lattice at B 7^ are shown in Table [2j 

It turns out that in the insulator phase, B < B c , the data for T](L) can 
very well be fitted by the exponentially decaying (in the L — > 00 limit) 
function 

rf\L) = Ce- L,L \ (18) 

where Lq and C are fitting parameters. The condensate in the infinite volume 
tends zero at B < B c , as expected. The fits are shown in Fig. [5] by the solid 
lines. The corresponding slopes are Lq = 0.42(2) fm and Lq = 0.90(8) fm for 
eB = and eB = 1.07 GeV 2 , respectively. 

At higher values of B, the condensate shows plateaux as L increases. 
We get the L — > 00 extrapolation for the condensate by averaging the data 
for i](L) at two largest values of the lattice size L (the horizontal dotted 
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lines in Fig. [5]). The extrapolated data - shown by the blue squares in 
Figure [3] - agrees quantitatively well with our small- volume analysis. The 
nonzero values of the extrapolated condensates are rj = 0.0046(4) GeV 2 and 
r) = 0.0095(1) GeV 2 for eB = 1.28 GeV 2 and eB = 2.14 GeV 2 , respectively. 
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Figure 5: The superconducting condensate t](L) vs. the lattice size L at fixed values of 
the magnetic field B. The dashed lines are shown to guide eye. The solid lines are the 



best fits of the data by the exponential function (181 



Thus, our numerical results support the theoretical prediction of Refs. 
|8] that the superconducting charged condensate of the p mesons forms spon- 
taneously at strong magnetic field. Using the simulations of the quenched 



QCD vacuum, we determined the critical magnetic field (16) which is re- 
markably close to the theoretical prediction ([I]). 

Finally, we would like to stress that theoretical calculations show that 
the condensate in the vacuum ground state should be an inhomogeneous 
function of spatial coordinates [TJ |8]. The ground state can be represented 
as a coherent static lattice of the topological (vortex-like) defects in the p- 
meson condensates, the p vortices, which are directed along the magnetic 
field axis. Qualitatively, the p-vortex state is very similar to the Abrikosov 
vortex lattices observed in the type-II superconductors in a background of a 
strong magnetic field [18] . 

It turns out, however, that in the QCD vacuum the energy gap between 
the lowest vortex energy state (given by a triangular vortex lattice) and 
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excited vortex lattice states is parametrically very small |T9], implying that 
the spatial lattice order of the vortex state may be destroyed by quantum 
(or thermal) fluctuations. The latter fact indicates that the actual vortex 
structures in the superconducting phase may resemble a much less ordered 
but persistent "spaghetti state", where the correlation functions, given by 
Eq. ([5]) and/or Eq. get additional suppression factors due to almost 
random vortex motion. The investigation of the detailed features of the 
superconducting ground state is currently underway [20J. 
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